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KECK/NIRC2 IMAGING OF THE WARPED, ASYMMETRIC DEBRIS DISK AROUND HD 32297 
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ABSTRACT 

We present Keck/NIRC2 K s band high-contrast coronagraphic imaging of the luminous debris disk 
around the nearby, young A star HD 32297 resolved at a projected separation of r = 0.3-2.5" (~ 35- 
280 AU). The disk is highly warped to the north and exhibits a complex, "wavy" surface brightness 
profile interior to r m 110 AU, where the peaks/plateaus in the profiles are shifted between the NE 
and SW disk lobes. The SW side of the disk is 50-100% brighter at r = 35-80 AU, and the location of 
its peak brightness roughly coincides with the disk's mm emission peak. Spectral energy distribution 
modeling suggests that HD 32297 has at least two dust populations that may originate from two 
separate belts likely at different locations, possibly at distances coinciding with the surface brightness 
peaks. A disk model for a single dust belt including a phase function with two components and a 5-10 
AU pericenter offset explains the disk's warped structure and reproduces some of the surface brightness 
profile's shape (e.g. the overall "wavy" profile, the SB peak/plateau shifts) but more poorly reproduces 
the disk's brightness asymmetry and the profile at wider separations (r > 110 AU). Although there 
may be alternate explanations, agreement between the SW disk brightness peak and disk's peak mm 
emission is consistent with an overdensity of very small, sub-blowout-sized dust and large, 0.1-1 mm- 
sized grains at ~ 45 AU tracing the same parent population of planetesimals. New near-IR and 
submm observations may be able to clarify whether even more complex grain scattering properties or 
dynamical sculpting by an unseen planet are required to explain HD 32297's disk structure. 
Subject headings: stars: early-type, planetary systems, stars: individual HD 32297 



1. INTRODUCTION 

Debris disks are signposts of planets and planet forma- 
tion (e.g. lWvattll200alKenvon and Broml ey 2008). Sup- 
porting this picture, the two stars with independently 
confirmed, direc tly imaged planetary syst e ms HR 8799 
and B Pictoris (iMarois et all [200l 120101 : ICurrie et all 
l2011at lLagrange et al.ll2010D are surro u nded by luminous 
debris disks dSmith and Terrild [l98l iRhee etail 120071: 
ISu et all I2009T K Similarly, Fomalhaut has a candidate 
planetary compani on located just in terior to the star's 
bright debris ring (IKalas et al.ll2008T ). 

In the absence of a directly imaged planet, re- 
solved imaging of debris disks may provide indirect ev- 
idence for a massive planet's existence, may help con- 
strain the unseen planet's properties, and thus can 
help iden tify promising targ ets for future direct imag- 
ing (e.g. iWyatt et all 119991 ) . For example, the in- 
clined ^)r_J|wjyTOed||coinrjonen^ debris 
disk (|Heap et alj I2000t iGoli mowski et all 120061) is likely 
due to the directly i maged planet (|Augereau et all 120011 : 
iDawson et all 120 111 ) and also provides an estimate for 
the planet's mass indepe ndent of planet cooling mod- 
els (jlagrange et all 120101 ). Dynamical sculpting by a 
planet/planets may explain the sharp inn er edge and 
pericenter offset of Fomalhaut 's debris ring (|Kalas et all 

1 NASA-Goddard Space Flight Center 

2 Steward Observatory, University of Arizona 

3 Space Telescope Science Institute 

4 NEXSCI, California Institute of Technology 

5 University of Wyoming 

6 Harvard-Smithsonian Center for Astrophysics 

7 Institute for Astronomy, University of Hawaii 

8 Keck Observatory 

9 MPIA-Heidelberg 



120051 IQuillen et all 120(35 IKalas et all 11)081) . Other de- 
bris disk structures may be due to non-planet processes, 
in particular interactions with the interstellar medium 
or perturbations from a nearby star, as has been pro- 
posed to explai n images of disks around HP 15115 and 
HD 61005 (e.g. IKalas et al.ll2007l; IHines eHI\^mi\) 

The nearby (d = 112±i§ pc; Ivan Leeuwenl 120071 ) A5 
star HD 32297 is another example of a young star sur- 
rounded by a lu minous, spatially-re solved, debris disk. 
At 30 Myr old (|Kalas et all 120051) . it is roughly co- 
eval with HR 8799 and may probe debris disk evolution 
at a stage just after they a re most collisi onally active 
(|Kenvon and Bromlevl 120081 : ICurrie et all 120081 120091) . 
like B Pic and HR 8799, HD 32297 has a large infrared 
(IR) excess emission due to circumstellar dust fi rst iden- 
tified from IRAS data. iSchneider et all (|2005| ) thus se- 
lected HD 32297 for Hubble Space Telescope(HST) NIC- 
MOS (F110W) coronagraphic imaging and resolved the 
disk out to an angular distance (from the star) of ~ 3.3" 
(~ 400 AU). HP 322 97 was subsequently re solved in the 
optical (IKalas! 1 20051 near-IR (1.6-2.2 fim. iDebes et all 
2009t iMawet et al.l 120091) . thermal infrared ( 10-20 p,m, 
Moerchen et all 120071: iFitzgerald etall 12007D . and mil- 



limeter fl.3mm lManess et alll2008l) . 

Previous work has claimed that HD 32297's disk struc- 
ture is sha ped by plane t sculp ting as well as non-planet 
processes. IDebes et all (|2009l ) identified an asymmetry 
in the disk scattering efficienc y between t he northeast 
and southwest sides (see also [Kalas 2005). They ar- 
gued that ISM sculpting explains this feature much like 
it explains s ome properties o f the HD 15115 and HD 
61005 disks. ISchneider et~aTI (f20Q5h identified a bright- 
ness asymmetry between the two disk sides, a feature 
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consistent with scu lpting by a massive planet (see also 
IManess et aHl2008ft . 

The two mechanisms, ISM sculpting and planets, are 
not mutually exclusive. New images of the HD 15115 
and HD 61005 disks reveal cleared inner regions and/or 
pericenter offsets, bot h of which are plausibly due to a 
plane tary companion (jRodigas et al.|[2012T: iBuenzli et al.l 
120101 ). Additionally, multiple debris belts, scaled ana- 
logues to the solar system's asteroid belt and Kuiper belt, 
are also likely planet signposts and m ay reside around 
HP 61005, HD 15115, and HD 32297 ([Fitzgerald et al.l 
12011 IManess et al.ll25ol iRodigas et al 1120121 ). 

To determine which mechanisms are responsible for 
shaping HD 32297's debris disk structure, we need new, 
high signal-to-noise images with which to derive p recise 
disk properties. Although Schneide r et al.l (|2005l ) iden- 
tify a disk brightness asymmetry consistent with planet 
sculpting, they caution that the disk brightness measure- 
ments close to the coronagraphic spot (r ~ 0.3-0.4") 
which provide the basis for this asymmetry are highly 
uncertain. PSF subtraction errors due to the completely 
opaque NICMOS coronagraphic spot may limit our abil- 
ity to conclusively identify disk structure at these small, 
speckle-dominated separations. Moreover, if the asym- 
metry identified a planet-induced density structure, it 
shoul d align with the mm emission peak (|Maness et al.l 
2008). However, it is not clear whether these asymme- 
tries are aligned and thus whether they identify small 
and large grains originating from th e same parent popu- 
lation. The Palomar/ifs image from Mawet et al. (2009) 
has limited spa tial resolution. While they did recover 
iSchneider et all s brightness asymmetry, higher spatial 
resolution observations could confirm and help clarify the 
physical origin of this and other features. For example, 
new data could identify breaks in the disk brightness pro- 
file that may reveal evidence for the multiple debris belts 
inferred from modeling unresolved IR data. 

To further clarify the nature of the HD 32297 debris 
disk, we present new coronagraphic imaging obtained at 
K s (~ 2 /im) with the Keck telescope on Mauna Kea, 
resolving the disk at angular separations of 0.3-2.5". §2 
describes our observations and extraction of the disk im- 
ages. In §3, we investigate basic disk properties (position 
angle, full-width half-maximum, and surface brightness) 
as a function of angular separation from the star. We 
then combine imaging with new, unresolved broadband 
photometry from the Spitzer Space Telescope and the 
WISE satellite to constrain the disk structure and iden- 
tify the location(s) of the disk emission (§4). Finally, 
we compare our analyses to those from previous work on 
HD 32297 (§5) and investigate the physical mechanisms 
responsible for sculpting the disk emission (§6). 

2. NIRC2 DATA 

2.1. Observations 

We imaged HD 32297 on UT November 20, 2011 with 
the NIRC2 camera mounted on Keck II using the K s fil- 
ter (A = 2.16 um) in the narrow camera mode (9.952 
mas /pixel lYelda et al.l[2010h with correlated double sam- 
pling. The Keck AO system delivered diffraction limited 
images with a FWHM of ~ 4.9 pixels (~ 49 mas). To en- 
hance our ability to extract the HD 32297 disk emission 
from the bright stellar halo, we centered the star behind 



the 0.6" diameter, partially-transmissive coronagraphic 
spot and used the "large hex" pupil plane mask. 

Our HD 32297 data consist of coadded 30s exposures 
with a cumulative integration time of 20 minutes and 
were taken through transit (HA = [-0.22,0.35]) in "ver- 
tical angle" or ang ular differential imaging mode (ADI, 
iMarois et all l20061 Over the course of these exposures, 
the parallactic angle changed by 36 degrees. While light 
cirrus caused fluctuations in the source transmission (as 
measured from the PSF core behind the coronagraphic 
mask) on the order of ~ 10%, observing conditions were 
otherwise stable. 

To flux calibrate the data, we measure the flux of the 
PSF core of HD 32297 as viewed through the partially 
transmissive coronagraphic mask and corrected for the 
extinction through the mask. To determine the extinc- 
tion, we use Keck/NIRC2 K s observations of HD 15115 
taken in August 2011 in photometric conditions with the 
same coronagraphic spot size that were flux-calibrated 
with a standard star HD 3029. From the August data, we 
measure the extinction through the 0.6" spot to be 6.91 ± 
0.15 mags, slightly less than but comparab le to extinction 
estim ates in K B for 1-2" spot diameters (Me tchev et al.l 
2009). Uncertainties in the aperture correction inside 
the coronagraphic spot due to scattered light from the 
spot edge limit the precision of our estimate. We did 
observe a photometric standard (p486r) 90 minutes after 
our HD 32297 imaging sequence. Although conditions 
became highly variable during the standard star obser- 
vations due to patchy cloud coverage, our absolute flux 
calibration implied from the least extincted standard star 
frames agrees with that derived from the coronagraph 
transmission to within 0.3 magnitudes. 

2.2. Image Processing 

Basic image processing and high-cont rast imaging 
techniques f ollow methods outlined in iCurrie et al.l 
(2010, 2011a, b). We employ standard dark subtraction 
and flatfielding corrections, identify and interpolate over 
hot /cold pixels , and apply t he di stortion correction de- 
termined from lYelda et al.l (|2010f ) . For image registra- 
tion, we exploit the fact that the PSF core is visible 
through the coronagraphic mask and is unsaturated. We 
register each image to subpixel (<7 con ~ 0.1 pixels) ac- 
curacy by determining the centroid position of the first 
image in the sequence and then determining the rela- 
tive offset between the reference image and subsequent 
images by solving for the peak in the cross-correlation 
function for each image pair. Finally, we subtract off a 
radial profile to remove the smooth seeing halo and make 
a second pass through the images to identify remaining 
bad pixels. 

To extract a detection of the HD 32297 disk, we adopt 
the Loc ally Optimized Combin ation of Images approach 
(LOCI ILafreniere et~aH 120071) . using an updated ver- 
sion of the LOCI-based code emplo yed in ICurrie et all 
(|2011bD . which will be detailed later (ICurrie et al.ll2012L 
T. Cu rrie, 2012 in prep.). Following iThalmann et all 
(|2011[) . we reduce the data using "conservative" LOCI 
settings more appropriate for extended sources (i.e. 
disks, not planets), where we consider rotation gaps S of 
1.5-5x FWHM, optimization areas (N A ) of 1000-3000x 
the FWHM area, optimization geometries g of 0.5-2, 
and subtraction annul (dr) of 5-10 pixels wide (see 
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lLafreniere et al.l 120071 for more details). To determine 
the signal-to-noise per resolution element of our disk de- 
tection, we convolve the image with a beam size equal to 
~ 1 FWHM and compare pixel counts to the standard 
deviation of counts within a 1 FWHM-wide annulusO. 
Finally, we correct for photometric / astrometric biases in- 
herent in LOCI-based processing by imputing fake disks 
into each registered image and comparing the predicted 
and measured disk properties (e.g. surface brightness, 
full-width half maximum, pos ition angle). Our me thod 
follows that first developed by iRodigas et al.l (|2012f ) and 
is described in full in the Appendix. 

2.3. HD 32297 Disk Image 

Figure [T] shows the image and signal-to-noise map for 
our "conservative" LOCI reduction assuming a rotation 
gap criterion of 5 > 2.5, Na — 3000, g — 1, and dr = 10, 
which balances our ability to detect the disk at small sep- 
arations by attenuating speckles (favoring smaller 5, Na) 
but without oversubtracting the disk (favoring larger <5, 
Na)- As evidenced by the signal-to-noise map, we detect 
the disk at SNR > 3 from r = 2.5" all the way to the 
edge of the coronagraphic spot at 0.3". The peak signal- 
to-noise per pixel is ~ 19. The disk emission along the 
midplane is more than 10cr significant at 0.85-1.4", while 
some regions of disk emission at 0.3-0.5" on the SW side 
are still more than 5-c significant. Visually inspecting 
the image and signal-to-noise map reveals significant disk 
structure. Most notably, the disk emission traces a dis- 
tinct "bow" shape, where the disk position angle clearly 
changes with radial separation. 

Furthermore, the disk exhibits a significant brightness 
asymmetry at small separations (r < 0.7"). Figure [2] 
redisplays our disk image with a different color stretch 
to better illustrate the brightness differences between the 
two sides of the disk. Interior to r = 0.35" (identified 
with a circle), the SW side is significantly brighter than 
both the NE side as well as any pixel value for residual 
speckles located at different azimuthal separations; the 
NE side has a peak brightness only slightly larger than 
the brightest speckle. Exterior to this separation, there 
are no residual speckles as bright as either side of the 
disk, and the SW is still clearly brighter at least out to 
r m 0.6" (yellow/red region on the SW side). 

Our Keck K s image agrees well with the previous 
best K s results, which were obtained with the extreme- 
A O Well- Correc t ed Su b-Aperture on Palomar presented 
by iMawet et al.l (|2009D using classical (not ADI) imag- 
ing. Convolved to the Palomar/WCS beamsize (Figure 
H bottom pa nel), our image strongly resembles that of 
Mawet et al.l 's. The disk appears highly asymmetric with 
the SW side being brighter extending all the way to the 
coronagraphic spot (0. 3" in our image com pared to 0.4" 
in theirs). As with the IMawet et alj (|2009l ) data, the NE 
side appears fainter and truncated. Moreover, the bright- 
est portion of the disk on the SW side roughly ov erlaps 
with the mm continuum peak (|Maness et~ayl2008h . 

3. ANALYSIS 



10 The focus of this paper is imaging and characterizing the HD 
32297 debris disk. While we do not present data reductions with 
methods more optimized for point source detection, we plan to do 
so in a future work 



Here, we investigate the HD 32297 disk geometry 
and surface brightness profile, following met hods simi- 
lar to those described in IRodigas et a l. (2012). We per- 
form analysis on the image shown in Figure [TJ The 
disk properties we report are corrected for photomet- 
ric/astrometric biases inherent in LOCI processing (see 
Appendix) . 

3.1. Disk FWHM 

To better assess the HD 32297 disk morphology, we 
measure the disk FWHM perpendicular to the disk's ma- 
jor axis as a function of stellocentric distance. First, we 
identify the brightest pixels at each radial separation for 
the NE and SW lobes, respectively. Next, we place a 5 
pixel by 21 pixel box centered on the brightest pixel and 
sum up the counts/pixel along each row of the box, pro- 
ducing a ID array of 21 values. Finally, we fit a Gaussian 
to this array, which yields the disk midplane location and 
the disk FWHM at that location. 

Figure [3] displays the disk FWHM as a function of stel- 
locentric distance for the NE (purple) and SW (green) 
sides. The errors correspond to the residuals of the Gaus- 
sian fits divided by the "throughput" for the disk FWHM 
as determined in the Appendix. On both sides, the disk 
FWHM steadily decreases from ~ 0.25" at r = 1.5" to 
w 0.15" at r = 1". Interior to r = 1", the disk FWHM 
fluctuates about a constant value, though this behavior 
is likely due to the difficulty of fitting a Gaussian pro- 
file in speckle-dominated regions, not bona fide clumpy 
structure. Beyond r = 1.5", where the disk emission ap- 
proaches the photon noise limit, the FWHM estimates 
fluctuate wildly. 

3.2. Disk Position Angle 

To quantify the "bow" structure easily seen in Fig- 
ure [1] and identify any additional warping, we calculate 
the disk position angle as a function of stellocentric dis- 
tance for both the NE and SW sides using the disk mid- 
plane pixel locations from the Gaussian fits described in 
§3.1. The position angle uncertainty at each radius re- 
sults from the difference between the Gaussian-fitted disk 
midplane location at the radius and a midplane location 
defined by the brightest pixel. Here we fo rmally assum e a 
systematic uncertainty of 0.009° from the lYelda et al.1 as- 
trometric calibration, although fitting errors always dom- 
inate. 

Beyond r = 0.9", both sides of the disk maintain a con- 
stant position angle, although they are misaligned by 3-4 
degrees. Between 0.3" and 0.9", though, the disk emis- 
sion on both sides curves towards the north, with the 
position angle decreasing on the NE side by more than 
20° and increasing by ~ 10° on the SW side. On the 
SW side, this curvature is not continuous, leveling off at 
232° at 0.5-0.6" before resuming at smaller separations. 
Beyond r = 1.6", the disk exhibits no obvious curva- 
ture/warping, though photon noise degrades the preci- 
sion of our estimates. 

3.3. Disk Surface Brightness Profile 

T o calculate the disk surface brightness (SB), we fol- 
low |Rodigai^F^L| (|2012l ) and determine the median sur- 
face brightness in mJy /arcsec 2 in a 18-pixel diameter 
circular aperture. We determine uncertainties in the me- 
dian disk surface brightness (also in mJy/ arcsec 2 ) at 
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a given angular separation in a way analogous to that 
which we used to determine the disk SNR/pixel. Specif- 
ically, at each pixel radius corresponding to a Gaussian- 
fitted disk midplane location we calculate the uncertainty 
in SB within non-overlapping circular apertures identi- 
cal in size to that which we use to determine the disk SB 
and covering all azimuthal angles exterior to the disk. 
We adopt the standard deviation of these SB estimates 
as the uncertainty in the disk SB at each radius. 

Figure |4] displays the surface brightness profiles for 
both sides of the disk. On both sides, the disk steadily 
brightens from r = 2.5" to r = 1". However, at smaller 
separations the profiles change dramatically. On the NE 
side, the disk has a roughly constant brightness at 0.5- 
0.9" before brightening from 3 to 6 mJy /arcsec 2 at 0.5" 
to 0.4". The SW side of the disk displays a similar be- 
havior: a nearly constant brightness at 0.7-0.9", a possi- 
ble slight dip in brightness at 0.7", and a sharp jump in 
brightness by a factor of 4 from 0.7" to 0.4". Thus, the 
disk SB profile appears "wavy" interior to r = 1". 

The two sides of the disk differ slightly in some other 
important respects. The locations of the peaks/plateaus 
closer to the star for the NE side than for the SW side. 
The shifted profiles are consistent with the dust ring be- 
ing located at different stellocentric distances between 
the NE and SW sides: a pericenter offset. The NE side 
also plateaus and may drop slightly at r ~ 0.3-0.35". 
Conversely, on the SW side the disk continues to brighten 
all the way to our inner working angle of 0.3". 

Moreover, our analysis confirms evidence for a bright- 
ness asymmetry between the NE and SW sides. These 
differences are significant even at small separations where 
we detect the disk at a lower SNR. Interior to 0.6", the 
median uncertainty in the disk SB is ~ 2.8 on the NE 
side and 4.6 on the SW side. While these uncertainties 
are large, they are significantly smaller than the bright- 
ness differences at r — 0.5-0.6", where the SW side is 2-3 
times brighter. At these angular separations, the +/- 1-cr 
range in SB between the two sides do not overlap. The 
SW side of the disk is 50% brighter at r - 0.3-0.4", al- 
though the large error bars for the NE side SB make this 
difference less statistically significant. 

We model the surface brightness profiles over r = 0.3- 
2.5" to derive power law indices a and b assuming a func- 
tional form of f(X) = aX b . As a first guess, we take 
the logarithm of the surface brightness, adopt uniform 
weighting to each datapoint, and estimate the power in- 
dices by fitting a straight line and deriving the slope. We 
then refine our power law index estimates by performing 
a Levenberg-Marquardt minimization, while considering 
measurement errors. To more precisely estimate the lo- 
cations of the power law breaks, we perform the above 
steps iteratively, varying the break locations and adopt- 
ing ones that minimize the reduced \ 2 between the data 
and the model. 

Table [T] summarizes our results. We fail to identify 
any power law to either lobe that quantitatively fits the 
entire radial extent of the disk. For the NE (SW) side, 
the goodness-of-fit value (R 2 ) declines to zero at r in ~ 
1.0" (0.9"). Assuming a single power law for both sides, 
we derive best-fit coefficients of a — 2.13, b — -6.01 at r = 
1.05-2.5" for the NE side and a = 1.72, b = -5.49 for the 
SW side at r = 0.95-2.5". However, the goodness-of-fit 



criterion for the NE side is low (R 2 ~ 0.75), indicating 
that at least this sid e may be best modeled as a broken 
powe r law (see also iSchneider et al.1 l2005t iDebes et al.l 
l200l . 

For the NE side, the best-fit indices assuming a broken 
power law are a = 2.23, b = -6.19 at r = 1.1-1.4" and 
a = 1.31, b = -5.13 at r = 1.4-2.5". For the SW side, 
the best-fit indices are similar: a — 1.78, b — -5.71 at 
r = 1-1.6" and a = 1.66, b = -5.33 at r = 1.65-2.5'EE 
Formally, our values for b have a large uncertainty since 
the goodness-of-fit criterion remains above 0.95 for a ± 
20% and b ± 0.5. 

We can compare our SB profiles to th ose from 
ISchneider etaLl (|2005l) and IDebes et~aTl (|2009j ) to under- 
stand how the profiles change with wavelength. In gen- 
eral, the shape of the SB profiles show good agreement 
with t hose derived from 1.1 fim data bv ISchneider et aLl 
(2005). Although they do not draw attention to any 
steep increase in SB at r = 0.5-0.7", their Figure 2 pro- 
vides some evidence for this feature, at least on the SW 
side. While we do not find evidence for a steep drop in 
SB at r = 0.3-0.5" for the NE side as they do, the SB 
does platea u at r = 0.4" and dr op slightly. The 2.05 [im 
profile from IDebes et al.l (|2009f ) does not extend to r < 
0.5-0.6", so we cannot know whether they too find evi- 
dence for a jump in SB at small separati ons. Our profiles 
are sig nific antly steeper t han th ose fromfSclmcidcr ~et al.l 
(|2005fl and IDebes et all (|2009f ) at r > 1 " obtained at 
shorter wavelengths (Figure [5|). 

4. DEBRIS DISK MODELING 

4.1. Scattered Light Modeling 

To understand the disk's grain properties and morphol- 
ogy, we compare our image to synethetic resolved disk 
models with a range of scattering properties. We mod- 
eled a number density distribution of dust for the disk in 
the following cylindrical form, which allows for a variety 
of disk morphologies: 

A(r,z)=exp ( ~ljjr j ex P {^^) ,r < rbrcak W 
( r 

= ) , r > r brC ak (2) 

\ '"break / 

where r Q can be interpreted as the location of a birth ring 
of planetesimals generating a collisional cascade of dust 
with a characteristic width a r . We assume a Gaussian 
scale height to the disk at all radii, and we allow for a 
power-law decrease in dust density beyond some radius 
'"break, thus allowing for combinations of ring-like and 
disk-like structures. We place dust particles at various 
distances from the star following the density distribution 
above and distribute them uniformly in the azimuthal 
angle (9). To derive the integrated surface brightness, 
we project the density distribution on a two dimensional 
plane after transforming the coordinates of dust particles 
to account for inclination and position angle and scale 
the result to the observed surface brightness of the disk 

11 While the inner and outer fitted radii reported here are the 
ones for which a power law fit is most applicable, we obtain similar 
results for slightly different radii: e.g. for Ti n = 0.95" and r ou t 
=1.65". 
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in the Keck images. The surface brightness, F, of each 
particle is determined from its distance from the star 
and its scattering angle ui assuming a two component 
Hcnyey-Greenstein phase function: 



F oc 
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A multi-component phase function for circumstellar 
dust may be favored. For example, the phase function 
of zodiacal dust in the Solar System has been modeled 
with multiple Henyey-Greenstein compone nts, includ- 
ing a significant backscattering component (|Hong et al.l 
I1985D . Recently, observations of HR 4796A's disk sur- 
face brightness as a function of scattering angle showed 
an incredibly flat inferr ed phase function for it s dust at 
scattering angles > 50° (D alle Ore et al. Il2011[ ) , and ob- 
servations of the protoplanetary disk HD 100546 at mul- 
tiple wavelengths require forward scattering grains where 
the ph ase function become s flatter at scattering angles 
> 38° (Mulder s et al.ll2012T) . 

Figure [5] displays the model that best reproduces the 
overall disk morphology and two rejected models. We 
attempted several possible structures and combinations 
of parameters, obtaining a good fit to the data with the 
model presented in Table [2] and shown in the left panel. 
A more rigorous exploration of parameter space and their 
possible degeneracies is beyond the scope of this paper. 

To account for its warped, bow-shaped appearance, the 
disk must contain some highly forward scattering grains 
at r ~ 110 AU, which cause a brightness asymmetry be- 
tween the front and back sides of the disk. At low scat- 
tering angles (small projected separations) the brightness 
asymmetry is more pronounced than at larger projected 
separations, causing a change in the midplane position 
angle. This t ype of "warping" structure is also seen 
in HD 15115 ()Debes et al J [20081: iRodigas et all 120121 ) . 
However, the sharp breaks in the surface brightness at 
~1" are hard to reproduce with typical single component 
Henyey-Greenstein phase functions (middle panel), re- 
quiring a flatter phase function at larger scattering angles 
(left panel). The surface brightness breaks could con- 
ceivably be reproduced by a superposition of two rings 
of isotropically scattering dust (right panel), but such 
a configuration does not give rise to the disk's warped 
appearance. We cannot exclude the presence of a sec- 
ond inner ring (i.e. at r ~ 35-50 AU) for our modeled 
dust phase function, especially if the dust in the disk at 
wider separations is slightly less forward scattering than 
we have modeled. This would allow the inner disk to 
dominate the surface brightness at small projected sepa- 
rations. 

Although our two-component forward scattering model 
is simple, it reproduces some key disk features. Figure [7] 
and [5] (left panel) compares this model to the observed 
SB profile. The model clearly succeeds in reproducing 
the disk warping at r < 1" and the surface brightness 
profile break/ turnover at ~ 110 AU. The model grains' 
strong forward scattering at small angles causes the disk 
to appear very bright again at small projected separa- 
tions (r < 0.6"). This feature agrees with our measured 
SB profile (Figure [8]), though taken at face value the 
model predicts that this brightness accelerates, whereas 



the measured profile appears to flatten at r ~ 0.4", espe- 
cially for the NE side. Furthermore, our scattered light 
model predicts that the disk emission should be roughly 
axisymmetric. However, at r < 0.6" the disk's SW side is 
significantly brighter and the profile breaks appear offset. 
The disk's asymmetric SB profile breaks could indi- 
~ fecate an , symmetry in the disk's distance from the star 
2g 2 <S8s W fufacXton of position angle: e.g. a pericenter offset. 
To investigate whether a pericenter offset better repro- 
duces the SB profile, we reconstruct a scattered light 
model with identical grain scattering properties as our 
two-component model but make the SW side of the disk 
5-10 AU closer to the star than the NE side (Figure |5J). 
This model predicts that the NE side's SB profile break 
starts at wider separations (r ~ 1.1" vs. 0.9" for the 
SW side) and that the SW side should be brighter at 
0.5-0.9", in agreement with the observed SB profiles. 

The pericenter offset model's fidelity isn't perfect: the 
SW side is significantly underluminous and the exact 
locations of the breaks are not well reproduced. How- 
ever, its success in reproducing the asymmetric SB pro- 
file breaks motivates more detailed scattered light disk 
modeling. In particular, all single ring models (regard- 
less of any pericenter offset) predict that the disk's sides 
should be equally luminous at r ~ 0.3-0.4". Our data 
indicate otherwise (the SW side is still brighter by ~ 
50%), suggesting the need for further modifications for 
our scattered light model. New, higher SNR images of 
HD 32297's disk will clarify the angular separation range 
over which the SW side is brighter and thus provide im- 
portant constraints for future disk modeling. 

4.2. Spectral Energy Distribution Modeling 

For a separate but related probe of the HD 32297 disk 
properties, we modeled the disk spectral energy distri- 
bution (SED) from point-source photometry including 
newly-available data from the Spitzer Sp ace Telescope, 
WISE mission, and AKARI satellite (iWerner et al.1 
l200l iWright et alJl20ia iMurakami et al.l 12007ft . Table 

[3] lists the photometric data we consider 

Our methods follow those outlined in iPlavchan et al.l 
(2009), where we identify the best-fit dust tempera- 
ture^), grain properties, and effective emitting areas in- 
corporating a downhill s i mplex algorithm ( "Amoeba" ) as 
described in lPress et al.l (|1992l ). We consider cases where 
the dust radiates like a blackbody such that r<iust,AU — 
(280K/T) 2 x^J L st ar/ Lq and where the dust's emissiv- 
ity scales with the effective grain siz e beyond a criti- 
cal wavelength (|Backman et al.l Il992f) . e oc A - ^, which 
can place grains of a given temperature at larger dis- 
tances. We solve for the disk model parameters that 
minimize the fit residuals relative to the flux (rms re i = 
J (X)"=o Ai/Fluxi) 2 /N). In all cases, we assume that 
the grain populations have a characteristic size: for sim- 
plicity, we do not consider a grain size distribution. 

Even though we include photometry not previously 
modeled, our fits are likely to be highly de generate (e.g. 
iManess et all l2008t iFitzgerald et al.1 120071) . Therefore, 
instead of identifying the single best-fit model and 95% 
confidence interval, we use several separate SED model 
runs to explore more focused questions about the HD 
32297 disk properties: 

• Model 1 — A Single Dust Population - Here we 
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assume that only one dust population contributes 
to the disk's IR-to-mm emission and allow /3 to 
vary. Our goal with this fit is to determine whether 
HD 32297's disk must have more than one dust 
population. 

• Model 2 Two Dust Populations, One Dust 
Belt at a Fixed Location - Here we include two 
dust populations and allow their emissivity laws to 
vary but fix them to the same location. We place 
the belt at r = 85 AU, or roughly interior to where 
the disk surface brightness profile begins to flatten. 
Our goal here is to determine whether more than 
one dust location is required. 

• Models 3 5 Two Dust Populations/Belts, 
Fixed Emissivity - Here we identify the best- 
fitting model with two dust populations in two dis- 
tinct belts, assuming either that the grains behave 
as blackbodies or have a A" 1 emissivity law. This 
approach follows that of iManess et al. (2008) and 
will allow us to assess whether their formalism pro- 
vides a better match to the photometry than as- 
suming one dust population/belt. 

• Model 6 Two Dust Populations/Belts, 
Variable Emissivity, Fixed Location - Here, 
we fix the locations of the dust to the peak of the 
surface brightness profiles: ~ 0.4" and ~ 1" or 45 
AU and 110 AU, which roughly cover the locations 
of the SB plateaus in the NE and SW sides. Our 
goal here is to see whether we can identify a good- 
fitting disk model that identifies the local maxima 
in the surface brightness profile as the locations of 
two separate debris belts. 

• Models 7-9 Three Dust Populations/Belts, 
Varia ble Emissivity, V ariable /Fixed Loca- 
tion - IManess et alJ ([2008?) suggest that three sep- 
arate dust populations are needed to fit the HD 
32297 SED. In one case, we fix all the belts to be 
at 85 AU and incorporate a third, warmer dust 
population varying in grain size and emissivity. In 
two cases, we fix the outer two dust belts at 45 AU 
and 110 AU, while varying grain size and emissiv- 
ity and incorporate a third warmer dust population 
varying in in radius, grain size, and emissivity. 

Figure [9] displays some of our SED modeling results, 
which are reported in Table [4] For all best-fit models, 
the stellar effective temperature is T* ~ 7870-7890 K 
and the star has little extinction (Ay < 0.03). Regard- 
less of our assumed particle emissivity law, a single dust 
population model (Model 1, top- left pa nel) poorly repro- 
duces the SED, as was found previously ([Fitzgerald et alJ 
l2007h . significantly underpredicting the flux density at 8- 
22 \xm. Formally, having two or three dust populations 
in one belt at 85 AU (i.e. top-right panel) significantly 
improves the fit (rms re i — 0.05-0.08). However, the re- 
quired grain sizes are too small (~ 4 nm-0.2 jim) to be 
realistic. Thus, the dust likely arises from more than one 
single-temperature dust population with given grain size 
and emissivity power law. 

Fits assuming blackbody grains but incorporating two 
dust components at different locations significantly im- 
prove the fit (rmsrei ~ 0.16; top-right) compared to 



a single dust component model. The fits imply that 
the dust responsible for 10-100 jim emission is at 1.2 - 
22 AU or r ~ 01-0 .2", b ut iFitzgerald et all rf200l . 
iMoerchen et al.l (|2007h . and IManess et all (|2008| ) show 
that the disk emission instead originates on scales more 
comparable to w 30-100 AU. The broad range of or- 
bital radii is indicative of the significant model degen- 
eracies involved. Assuming that at least one of the dust 
components has an emissivity of j3 = 1 yields fits with 
larger grains (0.3-3 /im) but worsens the fit (rms re i ~ 
0.17-0.21). Thus, if the dust disk geometries approxi- 
mate thin, isothermal rings and consistent of grains dom- 
inated by a single, characteristic size, the grains likely 
have emissivity power laws somewhere between and 1. 

The bottom panels of Figure [9] show that it might be 
possible for the debris rings to reside at locations equal to 
the scattered light brightness peaks (r ks 45 AU, 110 AU) 
and reproduce the IR to mm disk emission. A two dust 
population model (Model 6; bottom-left) reproduces the 
SED well at 1-8 /im and 17-1300 /im but underpredicts 
the disk emission at ~ 12 /im by ~ 50%. Starting with 
the Model 6 results and incorporating a third dust pop- 
ulation, we reproduce the SED if the dust is at ~ 14 
AU (171 K), although the emissivity power law for this 
population is unphysical (ft ~ 7.3)0 and the grain sizes 
for the coldest dust component are too small (Model 8). 
However, allowing the emissivity law and grain sizes for 
all three components to vary yields an excellent fit to 
the data (RMS ~ 0.06). The emissivity laws and grain 
sizes for this model (Model 9; bottom-right) are reason- 
able = 0.37-0.77; a = 0.20-0.70 /zm). Here, the third 
component is at ~ 1 AU with a temperature of w 430 
K and is responsible for the weak 8 [im excess and sub- 
stantially contributes to the disk emission at 12-17 \im. 
Thus, although it is not necessary, it is at least possible 
to identify the surface brightness profile peaks at 45 AU 
and 110 AU with separate thin dust rings responsible for 
broadband emission at 8-1300 /im provided that there 
exists an unresolved, warmer dust component located in- 
terior to 45 AU. 

We emphasize that the model fits we report in Table 
[4] are but some of many possibile fits, the range of which 
is further limited by our input assumptions about the 
disk. For example, because our model is set up only 
to consider infinitely thin, isothermal rings, it is unclear 
whether the warmer dust emission we identify originates 
in a separate ring or rather an annulus, the outer edge 
of which we see at wider separations. Given the ex- 
treme number of model degeneracies, the most we can 
say is that 1) the debris emission must originate from 
more than one population at 2) multiple locations and 3) 
among the many possibile configurations, the dust pop- 
ulations may coincide with the surface brightness peaks. 
To make further progress, we need high spatial resolu- 
tion, high SNR imaging of the disk at a wider range of 
wavelengths to derive much more stringent constraints on 
the dust sizes/scatterin g properties and location(s) of the 
dust population(s) (e.g. lDebes et al.ll"2009t iRodigas et al.1 
[2011) . 

12 This model run finds /3 ~ 7.3 as the "best-fit" model because 
the Amoeba code treats as a true free parameter, regardless of 
whether the value is physical, and exploits a gap in SED coverage 
at ~ fO /im to achieve a better fit. 
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5. COMPARISONS TO PREVIOUS HD 32297 RESULTS 

From optical coronagraphic imaging, IKalas et al.l 
(|2005T l find evidence for warping in the HD 32297's disk, 
especially for the NE side, where the disk at r > 500-600 
AU appears swept back compared to midplane regions 
at smaller angular separations. They attribute the warp 
to ISM sculpting of the disk as it moves south through 
ISM material. The warped structure we resolve curves 
in the opposite direction from that expected due to ISM 
sculpting (see im ages for HD 32297/15115/61005 from 
iDebes et al.ll2009D and does so for both sides of the disk. 
Furthermore, ISM sculpting should be prominent only at 
larger separations in our images, where the disk is weakly 
detected (r > 1.5"), because the grai ns need time to be - 
come entrained in the ISM flow (e.g. IDebes et al.l l2009). 
However, optical and near-IR imaging probe two very dif- 
ferent characteristic grain populations, and the size scales 
over which we resolve the disk do not overlap. Thus, our 
results are not in conflict. Rather, ISM sculpting may ex- 
plain the optical image but does not explain the near-IR 
image. 

We confirm and cl arify three major claims from previ- 
ous near-IR imaging (ISchneider et aL 2005; M awet et al.l 
I2009t IDebes et al. l l2009D. First, we ve rify t he brightness 
asymm etry found by ISchneider et all and iMawet et all 
(2009) at small (r ~ 45-55 AU) separations and show 
that it persists down to r ~ 35 AU, albeit at a 
lower statistical significance. Second, we recover a 
br eak / deviation in the SB power law at r ~ 100 AU found 
bv IDebes etall |g00g ). Our new results show that the 
SB profile at smaller separations does not follow an ex- 
tension of the power law describing the disk at larger sep- 
aratio ns: instead, it appears wavy. Finally, IMawet et all 
(2009]) claim that the NE side of the disk emission is 
truncated at r ~ 0.6". We find that the NE SB profile 
flattens (not brightens) from r ~ 1" to r ~ 0.5". As- 
suming that their coronagraph attenuates some flux at 
separations slighter greater than their inner working an- 
gle (0.4"), most of the NE side emission will be hidden 
from their view. In contrast, the SW side clearly begins 
to brighten at wider s eparations (r ~ 0.5-0.6") visible 
by IMawet et all ([20091 ) . Thus^ our results agree. 

Thermal IR im aging from iFitzgerald et all (|2007l ) and 
iMoerchen et al.l (|2007l ) shows that the HD 32297 disk 
has an inner clearing devoid of grai ns slightly larger than 
those we can probe with our data. IFitzgerald et al.l find 
that the disk exhibits a bilobed structure and most of 
its 11.2 [xm emission originates from r ~ 0.5-0.6", char- 
acteristic of a d ust ring with an inner edge at rs 60 AU. 
IMoerchen et al.l resolve the disk at 12 /mi and 18 /im out 
to slightly wider separations (r < 1.3"). Based on the 
disk's brightness temperature, they likewise find evidence 
for an inner clearing, albeit one that is slightly larger (r ~ 
80 AU) . The plateaus in our SB profiles co uld identify the 
inner boundaries of dust populations (e.g. iRodigas et all 
l2012f ) . Under this interpretation, our image is consistent 
with a dust belt truncated at r ~ 1" (w 110 AU) and 
a second one at w 0.4-0.6" (45-70 AU). Because these 
authors did not analyze their data in the same way it 
is difficult to compare their results between each other 
and their combined results against ours. Still, their re- 
sults and ours are qualitatively consistent with the HD 
32297 disk having at least one dust population truncated 



at small separations. 

In agreement with iManess et all (j2008[) 's mm study, 
we find evidence for at least two dust grain popula- 
tions, responsible for the mid-IR excess e mission and far- 
IR/su bmm excess, respectively (see also IMoerchen et all 
120071 ) . They are able to fit the SED from 1 to 100 \xm 
but not 1300 /im, a discrepancy they explain by adding 
a third, cold dust population. Our modeling generally 
fits the SED from 17 fim to 1300 (xm but slightly under- 
predicts emission at 12 /im, a discrepancy we can solve if 
there exists an additional, unresolved warm dust popu- 
lation. Our different results are likely a byproduct of our 
modeling assumptions: they allow the dust populations 
to be spatially extended but fix the particle emissivity 
law, whereas we assume the dust is confined in isother- 
mal rings but allow the emissivity law to vary. 

As noted in §3, t he peak of HD 32297's mm emission 
from IManess et all lies close to the bright inner disk re- 
gion on the SW side. While formally the position angle 
(« 46°) of the mm emission is offset by rs 10 degrees, 
it tracks the disk's position angle at wider separations 
much better. The disk's warped appearance at r < 1" is 
due to its strong forward scattering at small angles. If 
the disk's grains instead isotropically scattered starlight, 
its K s emission would lie almost directly on top of the 
mm peak. Thus, we identify the mm and near-IR bright- 
ness asymmetry as originating from the same location. 
Assuming the grains responsible for both the mm and 
near-IR emission are likely the result of collisions, they 
may trace the same parent population of planetesimals. 

5.1. Summary and Future Work 

Using Keck/NIRC2 K s coronagraphic imaging, we re- 
solve the HD 32297 disk at a high signal-to-noise from 
r = 0.3" to r = 2.5". We determine basic disk prop- 
erties (SB, FWHM, position angle), compare our image 
to disk models with a range of (grain) scattering proper- 
ties, and model newly-available, broadband photometry 
to provide a complementary investigation of HD 32297's 
circumstellar environment. We obtain the following ma- 
jor results: 

• 1. We discover that HD 32297's debris disk ex- 
hibits a prominent warped or "bow" -shaped struc- 
ture interior to r ~ 1 " (~ 110 AU). 

• 2. Our new data clarifies the disk's surface bright- 
ness profile at small separations. We find that it 
has a "wavy" profile interior to r ~ 110 AU with a 
plateau extending to r ~ 0.5-0.7" (55-80 AU) be- 
fore the disk brightens by factors of 3-4 at smaller 
separations. 

• 3. The disk exhibits significant asymmetries be- 
tween the two sides. The SW side is brighter at 
r ~ 0.3-0.6" by 50-150%, with the most statisti- 
cal differences being at r ~ 0.5-0.6". These sepa- 
rations are roughly where the disk's mm emission 
peaks and are consist ent with previous results (e.g. 
ISchneider et alll2005| ). Our new analysis identifies 
new asymmetries, revealing that the locations of 
the peaks/plateaus in surface brightness are likely 
shifted between the two sides, consistent with non- 
azimuthally symmetric structure. 
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• 4. A disk model with a flat phase function and 
strongly forward scattering grains where the dust 
ring is centered on the star reproduces the "bow" 
structure, marginally reproduces the "wavy" SB 
profile and fails to reproduce the NE/SW asym- 
metries. A disk model with a 5-10 AU pericenter 
offset reproduces the asymmetric SB profile breaks 
interior to r < 1", although its brightness asym- 
metry is limited to 0.5-0.9" and its match to the 
SB at wider separations is far poorer. Thus, dust 
scattering plays a critical role in explaining key ob- 
served disk properties, but it is unclear whether it 
explains all of the disk's properties we identify. 

• 5. HD 32297 must be surrounded by more than 
one dust population likely arising from different lo- 
cations in the disk. Although these populations 
need not be identifiable from our image, we can fit 
the disk SED by placing dust populations at the 
locations of the SB peaks from our image provided 
that there exists additional warm dust that we can- 
not yet resolve. 

• 6. The disk's brightness peak at r 0.4" coincide s 
with the peak mm emission (|Maness "etHI [20081. 
If the grains responsible for both peaks are the re- 
sult of collisions, they may trace the same parent 
population of planetesimals. 

In summary, the HD 32297 di sk appearance is 
broadly shaped b y ISM interactions (jKalas et al.l 120051 
iDebes et al.ll2009() at wide separations and by its grain 
scattering properties at small separations (this work). 
However, it is unclear whether either of these features 
by themselves explain the disk's SB profile and (espe- 
cially) the disk's brightness asymmetry. Furthermore, 
SED modeling provides evidence for multiple dust popu- 
lations, possibly multiple belts. To explain the disk's SB 
profile/multiple dust populations and (especially) bright- 
ness asymmetry, we may need to invoke other mecha- 
nisms. 

Planets can sculpt dus t into debris rings (e.g. IKalasI 
[2001 IQuillen et all [2001 IKalas et al.l 12008 ). Further- 
more, planets can trap dust into resonant structures, 
which can appear as bright, overdense regions like the 
bright ness peak imaged he r e and in the mm at r ~ 0.4" 
(e.g. iLiou and Zoold 119991 : iKuchner and Holmanl 120031: 
iStark and Kuchnerl 1200811 . As argued bv iWvatt et all 
([20061 ) and lManess et al.l (|2008| ). the detectability of the 
resonant structure may be wavelength dependent. Large 
grains producing mm emission are fragments of the col- 
liding planetesimals in resonance, but cannot be rapidly 
removed by radiative forces. Thus they can trace the par- 
ent body resonant structure. Small, (sub)-micron sized 
grains likewise may trace resonant structure since they 
are preferentially produced in the most collisionally ac- 
tive, highest density regions and are otherwise quickly 
removed by radiation pressure. Grains with intermedi- 



ate sizes producing emission at intermediate wavelengths 
(e.g. thermal IR) can be pushed out of resonance by ra- 
diation pressure/PR drag but are not small enough to 
be rapidly removed from the system (Wyatt et al. 2006; 
though see Kuchner and Stark 2010). Thus, it is possible 
that our near-IR image, when combined with the mm 
image, identifies planet-induced structure. 

Further near-IR and (sub)mm imaging is required 
to verify whether the SB profiles and brightness 
asymmetries are bona fide evidence for an embedded 
planet. The current state-of-the-art near-to-mid IR 
high-contrast imaging facilities like the Large Binocu- 
lar Telescope have already shown great promise for re- 
solving disk's like HD 32297's in scatt ered light and 
imag i ng self-luminous gas giant planets (Rodiga s et al.l 
IMa ISkemer et all [2012) and upcoming planet im- 
agers like SCExAO on Subaru, GPI on Gemini-South, 
and SPHERE on the VLT will be even more capa- 
ble (lEsposito et all 120 lit iMartinache and Guvonl 120091: 
iMacintosh et al.ll2008tlBeuzit et al.ll2008l 1. Moreover, the 
thermal IR is well suite d for imaging exop l anets with a 
wide range of ages (e .g. Marois et al. 201Q; ICurrie et al.l 
l2011bt iRodigas et al.l 120111 ). especially for stars like HD 
32297 whose bright, edge-on debris disk degrades planet 
sensitivity limits in the near-IR. Imaging with the Ata- 
cama Large Millimeter Array (ALMA) can resolve the 
HD 32297 debris disk up to factor of ~ 50 better than 
the CARMA observations reported bv lManess et al.l and 
thus will provide a far better probe of any planet-induced 
structure in the HD 32297 debris disk. 
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TABLE 1 

Surface Brightness Profile Power Law 
Indices 



Angular Separation (") a,b(NE) a,b(SW) 



Single Power Law 
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2.23,-6.19 
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1.31,-5.13 
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Note. — Surface brightness profiles for the 
NE and SW sides determined by a Levenberg- 
Marquardt least-squares fit to the data. We assume 
a simple power law functional form of f(X) = aX b . 
Power laws in agreement with the data at the 95% 
confidence limit include values for a that vary by ~ 
10% and values for b that vary by ~ 0.5 dex. 
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TABLE 3 
HD 32297 Photometric Data 



Filter /Wavelength Flux (mjy) crFlux (mjy) Source 



-O/0.4380 


Iood.74 


28.75 


1 YOH(J-ll(trans) 


V / 0.545 


on n ni 


24. 10 


1 Y L'rlU-ll(trans) 


J/1.235 


1341.80 


34.96 


zMAbo 


-ff/1.662 


913.48 


50.58 


2MASS 


A" s 2.159 


611.41 


11.95 


2MAbo 


3.37 




8.21 


WISE 


4.62 


151.28 


3.78 


WISE 


8 


65.18 


2.75 


IRSep 


11.2 


49.9 


2.1 


Fitzgerald et al. (2007) 


11.66 


53 


5.3 


Moerchen et al. (2007) 


12.08 


55.25 


1.20 


WISE 


16 


71.28 


2.53 


IRSep 


18.3 


90 


13.5 


Moerchen et al. (2007) 


22.19 


212.99 


5.55 


WISE 


24 


225.2 


4.82 


IRSep 


70 


850 


60 


IRSep 


90 


823.2 


116 


AKARI 


160 


<460 




IRSep 


1300 


5.1 


1.1 


CARMA/Maness et al. (2008) 



Note. — TYCHO-I1 (trans) refers to TYCHO-II catalog data transformed into 
the standard Johnson-Cousins photometric system. IRSep refers to the IRS en- 
hanced products dataset as queried from the NASA/IPAC Infrared Science Archive. 



TABLE 4 
HD 32297 SED Modeling Results 



Model ID RMS Rdust,i (AU) Rdust,2 Rdust,3 Pi P2 A «i (p m ) a 2 a 3 



1 


0.23 


85 






0.94 






0.095 






2 


0.08 


85 


85 




0.72 


0.39 




0.004 


0.197 




3 


0.16 


1.21 


21.66 












0.39 


0.41 




4 


0.17 


55.64 


527.56 




1 


1 




0.30 


0.99 




5 


0.21 


15.60 


97.50 







1 




0.91 


2.86 




6 


0.15 


45 


110 




0.80 


0.31 




0.39 


0.025 




7 


0.05 


85 


85 


85 


0.84 


0.61 


0.45 


0.006 


0.024 


0.725 


8 


0.09 


11 


45 


110 


7.27 


0.80 


0.31 


1.59 


0.39 


0.03 


9 


0.06 


1.1 


45 


110 


0.37 


0.77 


0.43 


0.70 


0.37 


0.20 



Note. — RMS refers to the fit residuals relative to the flux rms re [ = yj (5Z" =0 Ai/Fluxi) 2 /N, 
where N is the number of flux density measurements (18). Rduat refers to the dust ring's stellocentric 
dist ance, /3 is the particle emissivity power law, and a is the grain size in microns, where 27ra = A 
(sec Backinan et al. 1992). Values in bold are fixed for a given model run, whereas others arc 'fitted' 
values. 
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Fig. 1. — Reduced image (top) and signal-to-noise map (bottom) for our NIRC2 HD 32297 data. The colorbar depicting units for the 
image are in counts, whereas they range from to 9a for the signal-to-noise map. The central dark region identifies the coronagraphic spot 
(r = 0.3"). The panels have the same size scale. 
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Fig. 2. — (Top) Our image from the td]9 panel of Figure 3D redisplayed wittSQl different color s^Qetch to better illustrate the significance 
of our disk detection at the smallest separations (r =0.3- 0.6") and the brightness asymmetry. (Bottom) Our image resampled to the 
same spatial resolution as the Palomar/WCS image from Mawct ct al. (2009). The green cross identifies the position (and positional 
uncertainties) of the peak brightness in the millimeter Qylancss et al.ll200Sl l. Both panels are displayed in units of counts. 
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R (arc-sec) R (arc-sec) 

Fig. 3. — Disk FWHM (left) and position angle vs. angular separation for the NE and SW sides of the disk. The disk narrows at smaller 
angular separations. The two sides of the disk are offset in position angle by ~ 3-4°; the disk curves towards the north starting at r = 
0.9". 




0.5 1.0 1.5 2.0 2.5 

R (arc-sec) 



Fig. 4. — Surface brightness profiles for the two sides of the disk. In agreement with previous work (e.g. Schneider ct al. 2005), the disk 
exhibits power law breaks at r = 1.5" and r = 1.1". We identify a strong jump in surface brightness starting at r Ri 0.5—0.7". 
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Fig. 5. — Comparisons between our best-fit power laws to the HP 32297 K s surfac e brightness profile (SW side) and fitted power laws 
for l.f-2.05 fim HD 32297 data from Schneider ct al. (2005) and Dcbcs ct al. (2009). Our fits are generally much steeper. We find that 
the surface brightness profile interior to r ~ 1" cannot be fit by a power law. 
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Fig. 6. — Scattered light models incorporating different grain properties/disk geometries: (left) a two-component Henyey-Greenstein 
model with strongly forward scattering grains at small scattering angles but weakly scattering ones at larger angles (larger projected 
separations), (middle) a simple forward scattering grain model, (right) two isotropically scattering dust rings. The lefthand model best 
reproduces the disk SB profiles. The units are in counts. 
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R (arc-sec) 

Fig. 7. — Comparisons betwe en the mo deled and measured disk position angles (right). The triangle identifies the location of the mm 
brightness peak from Mancss ct al. (2008). We use the two-component model with the dust ring centered on the star and with a 10 AU 
offset. 
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Fig. 8. — Comparisons between the modeled and measured surface brightness profiles for a dust ring centered on the star (left) and one 
with a 5-10 AU pericenter offset (right; the SW side is 5—10 AU closer). Both models reproduce the wavy SB profile. The offset makes the 
disk model qualitatively reproduce the differences in the observed SB peaks/plateaus at r < 1", although it degrades the model's fidelity 
at r > 1 ", especially on the SW side, where it is substantially underluminous. 
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Fig. 10. — Images used to model biasing from LOCI processing. (Left) Our final image with a fake disk added. (Right) Our final image 
where we add the fake disk to each registered image and process the set of images with our pipeline. In both cases, we rotate the image to 
the parallactic angle of the first image in the sequence (PA ~ -15.06°), not to true north as we do in Figurc[T] 
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Fig. 11. — (Top panels) Surface brightness (left) and disk FWHM (right) vs. angular separation for the fake disk. (Bottom panels) Ratio 
of the "observed" (after processing) and expected surface brightness (left) and disk FWHM (right) vs. angular separation. The dotted 
lines identify power-law fits to correct our SB and disk FWHM measurements for biasing. 




Fig. 12. — Difference in position angle for the fake disk before and after processing. Our processing does not bias the disk astrometry by 
more than 1-2° at any separation. 
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APPENDIX 

CORRECTING FOR DISK PHOTOMETRIC / ASTROMETRIC BIASES WITH LOCI 

The LOCI-based PSF subtr action approach can bias the photometry and a strometry of point sources like planets and 
extended structure like disks (jLafreniere et al.l l2007t iThalman n et al.l 1201 If ) ■ For point source companions, imputing 
fake companions into registered images at a range of angular separations, processing these images, and then comparing 
the output to input fluxes and centroid positions of the fake sources corrects for these biases. Correcting biases for 
disks not viewed perfectly edge-on/pencil-thin is more difficult since we do not know, a priori, the disk's true FWHM, 
our ignorance of which also limits our ability to 'debias' other disk properties whose determination depends on the disk 
FWHM (i.e. surface brightness). Here we describe our method for mitigating LOCI-based photometric/ astrometric 
biases, which largely follows that of iRodigas et al.l (|2012f) . 

We first construct model disks each with midplane brightnesses (and thus signal-to-noise ratios) about twice that of 
the processed, real disk image and surface brightness profiles comparable to the real disk profiles beyond r — 1" . We 
consider two FWHM values of 0.2" and 0.25" at 1", or ~ 20-40% larger than the values we get between r — 0.3 and 
1" prior to applying the corrections we derive in this section (Figure [TTJ top panels). For both runs, we imput the fake 
disk into registered images, rerun our processing pipeline and then compare the output and input surface brightness 
profiles, disk thicknesses, and disk position angles. 

As shown in Figure [Til (middle and bottom panels), our processing minimially biases the thinner model disk (FWHM 
= 0.2") exterior to r = 1", reducing its surface brightness and thickness by no more than 20%. Interior to r = 1", the 
surface brightness and thickness drop to no less than 60-75% of their original values. Biasing for the FWHM = 0.25" 
disk (not shown) is only slightly more severe, resulting in an additional ~ 5% (10%) drop at r > 1" (r < 1") in SB 
and FWHM. For both fake disks, LOCI minimally biases the disk position angle measurements (Figure [T2"|) . 

The "observed" FWHM for the model with the initially thinner disk (FWHM =0.2") is ~ 10% smaller than the 
"observed FWHM for the initially thicker disk (FWHM = 0.25") and more similar to that which we find for the real 
HD 32297 disk prior to applying any bias corrections. Thus, we derive throughput corrections for the disk surface 
brightness and FWHM for the thinner disk, fitting a simple, unweighted power law to data between r = 0.3" and r = 
1.5". 

The deviations in throughput for individual points vs. our throughput function are as large as ~ 10-20%. However, 
our modeling errors for the disk SB and FWHM are larger in regions where biasing from LOCI is important and at 
r > 1.5-2" where the disk is intrinsically much fainter. Therefore, we leave a more detailed, robust calibration of 
disk parameters from LOCI-processed images to a future work where a higher SNR disk detection at r < 1" will help 
improve our fitting precision. 



